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Abstract 

We consider mixed integer linear sets denned by two equations involving two integer variables 
and any number of non-negative continuous variables. The non-trivial valid inequalities of such sets 
can be classified into split, type 1, type 2, type 3, and quadrilateral inequalities. We use a strength 
measure of Goemans [11] to analyze the benefit from adding a non-split inequality on top of the split 
closure. Applying a probabilistic model, we show that the importance of a type 2 inequality decreases 
with decreasing lattice width, on average. Our results suggest that this is also true for type 3 and 
quadrilateral inequalities. 

1 Introduction 

Wc deal with the mixed integer set Pj := {(x, s) G Z 2 x R" : x = f + YP j= i r j Sj} with / e Q 2 \ Z 2 and 
r J e Q 2 \ {0} for all j — l,...,n. We refer to / as the root vertex and to the vectors r J as rays. The 
motivation for analyzing P] is that it can be obtained as a relaxation of a general mixed integer linear 
set. Therefore, valid inequalities for conv(Pf) give rise to cutting planes for the original mixed integer 
set. A thorough investigation of this model can be found in Andersen et al. QQI2]- 

We consider the set i?/(r 1 , . . . , r n ) := conv({,s e M™ : / + Y^=i r " ' s j 6 ^ 2 }) which is the projection 
of conv(Pf) onto the space of the s- variables. A closed convex set B C R 2 with non-empty interior is 
lattice-free if the interior of B is disjoint with Z 2 . Furthermore, B is said to be maximal lattice-free if it 
is not properly contained in another lattice-free closed convex set. Any lattice-free polyhedron SCR 2 
with / in its interior gives rise to a function ip B : R 2 H> R + which is the Minkowski functional of B — f. 
We recall that for a convex set K containing the origin in its interior, the Minkowski functional || • ||k 
of K is defined by ||r||x := min{A >0:r£ XK} where r G M 2 . We call £)? =1 ip B (r j )sj > 1 the cut 
associated with B. This inequality is valid for Rf(r, . . . , r n ). Conversely, any non-trivial valid inequality 
for . . . , r n ) can be obtained from a lattice-free polyhedron. 

Every lattice-free rational polyhedron B C R 2 with / £ int(-B) is contained in a maximal lattice-free 
rational polyhedron B. It follows that ip B (r J ) < ip B (r 7 ') for all j = 1, . . . , n and the cut associated with 
B is dominated by the cut associated with B. This implies that the benefit from adding the latter cut on 
top of the split closure is at least as high as the benefit from adding the former. In this paper, we focus 
exclusively on cuts associated with maximal lattice-free rational polyhedra. 

A classification of planar maximal lattice- free closed convex sets was given by Lovasz [14] . The refined 
classification which is stated below can be found in [TU] . 



Proposition 1.1 ( jlOl 114] ^1 Let B be a maximal lattice-free closed convex set with non-empty interior. 
Then B is a polyhedron and the relative interior of each facet of B contains at least one integer point. 
In particular, B is one of the following sets (see Fig. fJJ). 

/. A split {(xi, X2) S R 2 : b < a±Xi + 02X2 < b + 1} where ai and 02 are coprime integers and b is an 
integer. 

II. A triangle which in turn is either 
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(a) a type 1 triangle, i.e. a triangle with integer vertices and exactly one integer point in the relative 
interior of each edge, or 

(b) a type 2 triangle, i.e. a triangle with at least one fractional vertex v, exactly one integer point 
in the relative interior of the two edges incident to v and at least two integer points on the third 
edge, or 

(c) a type 3 triangle, i.e. a triangle with exactly three integer points on the boundary, one in the 
relative interior of each edge. 

III. A quadrilateral containing exactly one integer point in the relative interior of each of its edges. 




(a) Split (b) Type 1 triangle (c) Type 2 triangle (d) Type 3 triangle (e) Quadrilateral 



Figure 1: All types of maximal lattice- free convex sets in K 2 with non-empty interior. 

By Proposition ll.il it follows that i?/(r 1 , . . . , r n ) has three types of non-trivial valid inequalities: split, 
triangle, and quadrilateral inequalities named after the corresponding two-dimensional object from which 
the inequality can be derived as an intersection cut [5]- Triangle inequalities are further subdivided 
into type 1, type 2, and type 3 inequalities (see Fig. QJ. Note that a non-trivial valid inequality for 
Rf(r , . . . , r n ) can correspond to more than one maximal lattice- free polyhedron. 

We are interested in the quality of cuts associated with the maximal lattice-free polyhedra described 
in Proposition 11.11 We use the strength measure introduced by Goemans Ql] to evaluate the quality 
of a cut. Basu et al. [7J assess the strength of split, triangle, and quadrilateral inequalities in a non- 
probabilistic setting. They show that split and type 1 inequalities may produce an arbitrarily bad 
approximation of Rf^r 1 , . . . , r n ). On the other hand, type 2 or type 3 or quadrilateral inequalities deliver 
good approximations of Rf(r , ... , r n ) in terms of the strength. This, however comes with a price. Up to 
transformations preserving the integer lattice, there is only one split and only one triangle of type 1, but 
an infinite number of triangles of types 2 and 3 and quadrilaterals. Therefore, it can be expected that 
for real instances an approximation by adding all these cuts is hard to compute. From a more practical 
point of view, one is interested in approximations of the mixed integer hull that one can generate easily. 
Current state-of-the-art in computational integer programming is to experiment with split cuts and the 
split closure (see e.g. [SUB]). This is the point of departure of our theoretical study. 

The aim of this paper is to shed some light on the question which average improvement a non-split 
cut gives when added on top of the split closure. For that, we take any maximal lattice-free triangle or 
quadrilateral B, and we investigate all potential sets Rf{r l , . . . , r n ) such that we can generate a valid 
cut from B. For this, it is required that / is in the interior of B. We vary / uniformly at random in 
the interior of B. This is our probability distribution. For each particular B and / £ int(_B) we let n 
and r 1 , . . . , r" attain arbitrary values. We compute a lower bound on the probability that the strength 
of adding the cut associated with B on top of the split closure is less than or equal to an arbitrary value. 

As a conclusion from our probabilistic analysis we obtain that the addition of a single type 2, type 3, 
or quadrilateral inequality to the split closure becomes less likely to be benchcial the closer the lattice-free 
set looks like a split, i.e. the closer its geometric width is to one. 

We think that this complements nicely the analysis in [7J: they construct sequences of examples 
in which cuts from triangles of types 2 and 3, and quadrilaterals cannot be approximated to within a 
constant factor by the split closure. The approximation becomes worse as the triangles and quadrilaterals 
converge towards a split. From the results in our paper it follows that this geometrically counterintuitive 
situation occurs extremely rarely. 

2 Probabilistic model and main results 

Often, the quality of cuts is measured according to their worst case performance. In this paper we apply 
a stochastic approach by introducing a probability distribution on all possible instances. See [5] and [T2] 
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for other approaches on how to apply a probabilistic analysis to mixed integer linear sets. 

Our aim is to evaluate the benefit from adding a single cut associated with a maximal lattice-free 
rational triangle or quadrilateral on top of the split closure. For that reason we apply the following 
strength measure. 

A non-empty polyhedron P = {a; G R" : Ax > 6} is said to be of covering type if all entries in A and 
b are non-negative and ^ P. Let a G (0, +oo). We call aP := {x : ax G P} the dilation of P by the 
factor a and define aP := R" if a = +oo. Note that P C aP for a > 1. Let Q C R™ be any convex 
set such that P C Q. The strength of P with respect to Q, denoted by t(P,Q), is the minimum value of 
a > 1 such that Q C aP. 

Here, the underlying idea is that a high strength means high benefit. In other words, let P and Q be 
two relaxations for a mixed integer set M such that P G Q. The larger t(P, Q) the better is the gap to 
M closed by P compared to that of Q. 



Lemma 2.1 Let P := {x G R™ : ajx > bi for all i = 1, . . . , m} be a polyhedron of covering type 

and Q C R™ any convex set such that PCQ. Then, 

t(P, Q) = max ( . , t~~ 7^7 : 6 * > °1 ■ 

i=i,...,m ^mf{a/a; : x G J 

// inf{a^2; : x G Q} = /or some i G {1, . . . , m}, where bi > 0, then t(P, Q) is defined to be +oo. 

In the following we use BJl instead of i?/(r 1 , . . . , r n ) for simplicity. Note that RJ is a polyhedron |15j . 
Let be the set of all maximal lattice-free rational polyhedra in R 2 which contain / in the interior and 
let S (resp. 71, 72, 7-j, Q) be the subset of containing all splits (resp. triangles of type 1, type 2, type 3, 
quadrilaterals). Observe that R'J = {s G R™ : £™=i i> B (r j )sj > 1 for all B G 5 U 71 U % U 75 U Q}. For 
a non-empty set £ C f2 we define to be the intersection of all cuts associated with the polyhedra 

in C and the trivial inequalities Sj > for all j = 1, . . . ,n. In the remainder of the paper, C will always 
be iS or S U {B} for some B G Cl\S. In this case C(Kj) is a polyhedron [S] of covering type. 

In order to evaluate the gain from adding a triangle or quadrilateral inequality to the split closure we 
apply the following procedure. We fix a maximal lattice-free rational triangle or quadrilateral B G fl \5, 
but we allow / to vary in its interior. We will show that by considering a specific set of rays which 
depends only on /, we obtain an upper bound on t(J 7 (R , j),S(R''j)) where J = 5U {B}. Depending on 
where / is located this bound may differ. By varying / over the entire area of B we compute the area 
for which the bound is below a certain value, z say, and compare it to the area of B. This gives a ratio 
which is a lower bound on the probability that the strength is less than or equal to z. In turn, 1 minus 
this probability is an upper bound for the chance that B improves upon the split closure by a value of 
more than z with respect to the strength. 

Let B G O \ S and let T = S U {£>}. The following observation follows easily from Lemma \2. II 



Observation 2.1 

t(F(R]),S(R])) = ^ n{E n =i ^ B(rJ> .. sG5(fl „ )} . 

Since B, /, and r 1 , . . . ,r™ are rational, we can assume that the rays r 1 , . . . , r n are scaled such that 
the points / + r J are on the boundary of B (see [7] for an explanation why this assumption is feasible) . 
In the course of the paper we will deal with optimization problems of the following type: 

mmJ2^ jB ( rj ) s j s -t- seS{Rj) (1) 
j'=i 

for some B G £1 \ S. Scaling the rays as described above implies that the objective function becomes 

En 

We define a corner ray to be a ray r 3 where the point f + r J is a vertex of B. For now, assume that 
all the corner rays (and only those) are present. We will explain later why we can make this assumption. 
This implies that the corner rays are fixed once / and B are chosen. We assume a continuous uniform 
distribution on / in the interior of B. Given z G R, z > 1, we define P B (z) to be the probability that 
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t(J r (RJ),S(R 1 ^)) is less than or equal to z for / varying in the triangle or quadrilateral B, i.e. 

PB{Z) := A^{B) J mHR]),S(R]))<z}df, 

/Sint(S) 

where Area(i?) is the area of B and 1 is the indicator function. 

Let us now argue why we can assume that all corner rays of B are present in order to obtain the desired 
bound. Assume k > 1 corner rays r n+1 , . . . , r n+k are missing and consider R n ^ +k :~ conv({,s G : / + 

E"=i fc r3g j e Z 2 }). We now apply Observation O to infer t{F{R n f +k ) 1 S{R n f +k )) > t(T(R]), S(KJ)) <^=> 

min{X)™=i Sj : s € <S(i?^ +fc )} < min{^™ =1 Sj : s £ S(R'j)}. The latter inequality follows from the fact 
that an optimal solution s for the latter minimization problem implies a feasible solution for the former 
minimization problem by setting Sj = for all j = n + 1, . . . , n + k. This is because if there would be a 
split which cuts off (s, 0) from iS(i?^ +fc ), then the same split would cut off s from iS(i?^). It follows that 

l{t(F(R™ +k ),S(R 7 } +k )) <z}< l{t{F{R r }),S{R n f )) < z}. Thus, by adding the corner rays we obtain a 
lower bound on P B (z). 

We now explain why we can assume that only the corner rays of B are present in order to compute 
t(J r (R 7 ^ +k ) , S(R^ +k )) . In [7] it is shown that, if all corner rays are present and under the assumption of 
scaled rays, ([1} reduces to the problem where only the corner rays are given, i.e. the objective function 
is X)j=i s jj where {r 1 ,...,r fc }, k & {3,4} is exactly the set of corner rays for the given triangle or 
quadrilateral B. By Observation IO this implies t(J 7 (R , } +k ),S{R , } +k )) = t(J : (R k f ),S{R k f )). Thus, in the 
remainder of the paper, in order to obtain the desired bounds, we assume that the set of rays consists of 
exactly the corner rays of B. 

An important value in our analysis is the so-called lattice width. Let B C R 2 be a maximal lattice-free 
set with non-empty interior. The lattice width of B is defined to be 

w(B) := min <maxji r i-min« T i 

uGZ 2 \{0} [xeB xeB 

It is shown in |13| that 1 < w(B) < 1 + Furthermore, applying an affine unimodular transformation 

(i.e. a transformation preserving Z 2 ) to B changes neither the lattice width nor the strength. In the 
remainder of the paper, we informally call B flat whenever w(B) is sufficiently close to 1. Our main 
results arc summarized below. 




Theorem 2.1 (probabilistic strength - lattice width relation) LetTi be a triangle of type i e 
{1, 2} and w := w(T2). Then, for any z > 1, we have 



I. 



if 1< z < §, 

U(^r) 2 if§<2<2, 

if 2 < z < +oo. 




oo, 



With 01 - (z-™)(2v>z-v>-z) and ao - (w-lY(z-lY-l 
wlui yi — w 2 (z—l) 2 u,nu, </2 



2 (z-l) 2 i> z — w 2 (z-l) 2 

Theorem I2.1IIIII shows that for any given z > 1, P T2 (z) tends to 1 if wiT^) converges to 1, i.e. the 
probability that a fiat type 2 triangle improves upon the split closure by a value of more than z goes to 
0. This will be explained in further detail in Section [SJ 

The analysis of a type 3 triangle T3 and a quadrilateral Q turns out to be more complex. We did not 
succeed in putting P Ti (z) (resp. pQ{z)) into direct relation to w(Ts) (resp. w(Q)) and z only. Instead 
we parametrize T3 and Q in terms of the coordinates of their vertices. Using this more complicated 
parametrization we derive formulas for P T3 (z) and PQ(z). Then we discretize the coordinates of the 
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vertices and evaluate the formulas with respect to our discretization. This qualitatively leads to the same 
conclusion as before: if T3 and Q converge towards a split (meaning the lattice width converges to 1), 
the probability that the strength of the associated inequality is less than or equal to z tends to 1. We 
refer to Sections [5] and [7] for the corresponding formulas. 



3 Type 1 triangles 

By an affine unimodular transformation, we assume that the type 1 triangle T\ is given by T\ = 
conv{(0,0),(2,0),(0,2)}. Let R 1 := int(conv{(l, 0), (0, 1), (1, 1)}), R 2 := int(conv{(0, 0), (1, 0), (0, 1)}), 
R 3 := int(conv{(0, 1), (1, 1), (0, 2)}), and i? 4 := int(conv{(l, 0), (1, 1), (2, 0)}). 

Note that int(Ti)\u| =1 i?j is a set of area zero, so it can be neglected. For given f € R% UR2UR3UR4 
the three corner rays are r 1 = {-f lt -/ 2 ), r 2 = (2 - f u -/ 2 ), and r 3 = (-/1, 2 - / 2 ). Let J = 5U {T x }. 
In [7J it is shown that 

'2 if /e Ri, 



3-/1 



2-/i -/: 



if/eft, 



if / G R 4 . 



t(F(R}),S(R})) 

For the sake of a shorter notation we write l{t < z} instead of l{t(J-(R'j),S(R : j)) < z}. Then, 

1 

Area(Ti) 



We compute the four integrals separately. For that we need to check when the corresponding functions in 
Q attain a value which is less than or equal to z. Assume / G R±. Then ff eRl < z}df is if z < 2 

and is \ if z > 2. Assume / G R 2 . We have 3 ~^~^ < z <^> /1 + / 2 < and thus J /efl . 3 l{i < z}df is 
if z < |, is I if z > 2, and is \{ 2 ^Er) 2 otherwise (here we used the fact that < /1 + / 2 < 1 in R 2 ). 
Assume / G i? 3 . Then < z <^ / 2 > j^-. Hence, J /ei?3 l{t < is if z < |, is \ if z > 2, and 
is \{ 2 ^Zi ) 2 otherwise (here we used 1 < f 2 < 2). Finally, the case / G i?4 is analogous to the previous 
case. Since Area(Ti) = 2 it follows 

if l<z < §, 

^WHK^fr 1 ) 2 if I < ^ < 2, 

1 if 2 < z < +00. 



4 Proof strategy 

The analysis of the strength of type 1 triangles is quite easy since the split closure is known: it is always 
defined by a subset of the three split inequalities corresponding to the splits whose normal vectors are 
the normal vectors of the facets of the type 1 triangle (see [Jj for a proof). Note that w(Ti) = 2 and that 
this value is attained by precisely these vectors. 

Using the split closure for triangles of types 2 and 3, and quadrilaterals would result in too complicated 
formulas. Thus, we choose another strategy. Instead of using the entire split closure we will take only one 
well-chosen split inequality and therefore obtain lower bounds for the desired probabilities. Let B G fi\iS. 
The split inequality which we choose will depend on the location of / in the interior of B. For that we 
partition B into regions R± , . . . , R p and select a single split for each region. The basic idea is to choose a 
split that contains / in its interior and such that the normal vector of the split is a potential candidate for 
an integer vector for which w{B) is attained. Our candidate splits will always be among the vectors (1, 0), 
(0, 1), and (1, 1) since we apply an affine unimodular transformation to B to bring it into an appropriate 
form (see Sections [SJ H2 and [JJ for details). 

We now show that our simplification of using only one split inequality instead of the split closure 
leads to a lower bound for P B (z). Let T = S U {B} for some B £ !1\iS and let L G S be arbitrary. 
Then, by Observation [23 ^(i^), S(R k f )) = (min{^ =1 Sj : s G S(i^)})- 1 < (min{^ =1 Sj : Sj > 
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for j € {l,...,k} and E?=i ^V>i > l})" 1 = t({E?=i > l,Ej=i lW)*j > 1} H > 

for j e {1, . . . , fc}}, {Ej=i ^ L (r j )s 3 > 1} n {.Sj > for j € {1, ... , fc}}), since we assumed ^ s (r 3 ) = 1 
for all j = 1, ... , fc. 

For ease of notation we denote the latter by i(B, L), which is the strength of the polyhedron obtained 
by adding the cuts associated with B and L with respect to the polyhedron obtained by just adding the 

1%S{R)) 



cut associated with L. It follows l{t{F(R k f ), S(R k f )) < z} > l{t(B,L) < z} and therefore 



pB ^^A^B)t J l{mL Rj )<z}df, (4) 

where Lr, is the single split which is used in region Rj to approximate the split closure, for j = 1, . . . ,p. 
In the following, for simplicity, we write l{i < z} instead of l{i(B,L) < z} whenever B and L are clear 
from the context. In order to compute i(B,L) we need to solve an optimization problem of the type 

min si + • — h s& 

s.t. ip L (r 1 )s 1 + --- + Tjj L (r k )s k > 1, 
Sj > for j = 1, . . . , k. 

The values ^(r 1 ),.. .^ L {r k ) are called the coefficients of the split L. In general, for a split L = {i6 
R 2 : [tt t /J < ir T x < [~7r T /] }, 7r <E 1? \ {0}, containing / in its interior and with corresponding inequality 

Ej=x i jL ( r: ') s j > 1, it is easy to verify that 

if tt^V = 0, (5) 

1 ,^ rJ T , if ttV < 0. 

Thus, we just need the normal vector 7r of the split L to compute the split coefficients ^(r 1 ), . . . , ip L (r k ). 

In Scction[5]wc will explain our computations for triangles of type 2 in detail. Since the computations 
for quadrilaterals and triangles of type 3 in Sections [6] and [7] give no new insights we will only state 
intermediate results there. 



5 Type 2 triangles 

By an affine unimodular transformation, we assume that the type 2 triangle T2 has one facet containing 
the points (0,0) and (1,0), one facet containing (0,1), and one facet containing (1,1). Furthermore, 
one vertex a = (01,02) satisfies < a\ < 1 and 1 < a-i- Thus, the other vertices are (— a a ^_ 1 , 0) 
and ( ""^"j 1 , 0). We assume that a is arbitrary but fixed and treat it as a parameter in the subsequent 
computations. We decompose T 2 into six regions: R\ := int(conv{(0, 0), (01, 0), (01, 1), (0, 1)}), R2 := 
int(conv{( ai ,0),(l,0),(l,l),( ai ,l)}), R 3 := int(conv{(0, 0), (0, 1), (-^ T , 0)}), R 4 := int(conv{(l, 0), 

(1,1), (S L ,0)}), Rb := int(conv{(0,l), (01, 1), (01,02)}), and R e := int(conv{(oi, 1), (1, 1), (a x , a 2 )}) 
(see Fig.©. 
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Figure 2: Decomposition of a type 2 triangle. 

For given / <E L)j =1 Rj the three corner rays are r 1 = (— — /lj — A), ?* 2 = ( ~ Ai — /a)' an< ^ 
7' 3 = (ai — /i,a2 — A)- Furthermore, Arca(T2) = 2/n^-i) ' ^ ur assum pti° ns on ^2 imply that w(T2), as 



G 



denned in ©, is attained for u = (1,0) or u = (0, 1) (see @|). Thus, w := w(T 2 ) = min{ci2, ^zt) which 
implies w = a 2 if a 2 < 2 and to = a °^ 1 if 02 > 2. In either case we have Area(T2) = 2 (Z-i) • 

In regions R3 and R4 we use the split Si := {(xi,X2) £ M 2 : < X2 < 1} and in regions -R5 and i?6 
we use the split S2 ■= {(xi,X2) £ K 2 : < x\ < 1}. In regions R\ and R 2 we choose either Si or S 2 . We 
will use the one whose normal vector attains w, i.e. we choose Si if a 2 < 2 and 5*2 if a2 > 2. Note that 
our choice of the single split which is taken instead of the split closure is such that it attains minimal 
lattice width with respect to T 2 . Such a split has another nice property: among all splits that contain / 
in their interior we select one which covers most of the area of the triangle. 

5.1 Regions R 3 and i? 4 

Let / £ R3 U i?4. In order to compute i(T 2 , Si) we have to solve the optimization problem 

min Si + s 2 + S3 

s.t. ij Sl (r 1 )s 1 +iJ Sl (r 2 )s 2 + ^ Sl (r 3 )s 3 > 1, (6) 
Sj>0 for j = 1,2, 3. 

Since we use Si we have 7r = (0, 1) and since / £ R3UR4 it holds < ir T f < 1. Thus, using ([5]), we obtain 
^(r 1 ) = ip Sl (r 2 ) = 1 and -tjj Sl (r 3 ) = -fr^f ■ Therefore, an optimal solution to © is s* = (0,0, ^Er)r) 
with optimal objective value ^75- • It follows t{T 2 ,Si) = ^f^ 2 - for / £ i? 3 U i? 4 . 

5.2 Regions i? 5 and i? 6 

Let / £ R5. To compute t(T 2) S 2 ) in this case we solve the optimization problem min Y)j—-\ s j such 

that J2 3 j=ii' S2 ( rj ) s j > h s j > for 3 = M; 3 - Applying © we obtain ^(r 1 ) = h+ ]f* , ^(r 2 ) = 

"V-/! — , and if; S2 (r 3 ) = "^Zfl ■ Hence, the optimal solution is the minimum among {(V* v* 1 )) -1 ) 
(i/; S2 (r 2 )) -1 , (-0 s2 (r 3 )) -1 }. Using our assumptions on the variables and parameters it follows (^^(r 3 ))^ 1 
> 1 and (^(r 1 ))- 1 < 1 for i £ {1,2}. One easily verifies that (■0 s2 (r 1 ))" 1 < (ip S2 (r 2 ))" 1 01 > /1 

which is satisfied in R$ by assumption. Thus, t(T2, S 2 ) = — 2 ~ 1 for / £ i?g. 

a 2~°l _y 1 

Let / £ i?6- By symmetry, i.e. a± — > 1 — ai and /x — »• 1 — /1, we obtain t{T 2l S 2 ) = a2 ilj — for 
/ £ Re- 

5.3 Regions R\ and R 2 

First assume a 2 < 2 and use the split Si. Let / £ R\ U R 2 . The associated optimization problem 
is min Y?j=i s j s -t. Y J ] = i'4> 8l {f)sj > !> s i > for J = 1.2,3 with ^(r 1 ) = ip Sl {r 2 ) = 1 and 
ip Sl (r 3 ) = 2 fEjf- Hence, i(T 2 ,Si) = f or f eRiUR 2 and a 2 < 2. 

Now assume that 02 > 2. We use the split S2. Let / £ R\. The solution of the optimization problem 

min Ej=i s j s-t. T, 3 j=i^ S2 {r j )sj > 1, Sj > for j = 1, 2, 3 with ?/> S2 (r 1 ) = ^jf* ,ip s *(r 2 ) = "\l^ h , 
and -0 S2 (r 3 ) = ^f^- is the minimum of the set {(ip S2 (r 1 )) -1 , (V ,S2 (^ 2 )) -1 } as (ip S2 (r 3 ))^ 1 > 1. It is easy 

to verify that (i/> S2 (r 1 ))" 1 < (V> S2 (r 2 ))" 1 ^ ai > / x . Therefore, we have i(T 2 , S 2 ) = h+ jf^ for / £ R 1 
and a 2 > 2. 

Finally, assume / £ R 2 . By symmetry, we obtain t(T 2 , S 2 ) = °' 2 1 lf l — for f <E R 2 and a 2 > 2. 
The following table summarizes the function i(T 2 ,Si), for the corresponding i £ {1,2}. 
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i{T a ,Si) 
for a,2 < 2 



Location of / 




5.4 Approximation for P T2 (; 



We compute the integrals Jf eR . l{i(T"2, 5;) < z}df for j G {1, . . . , 6} and the corresponding split Si or 
5 2 which we used above. For simplicity, ff eR . '■= ff e n- ^ z l c ?f f° r J = 1, ... , 6. 

Let / G i? 3 . Then 1{F < z] = 1 <S> < z f 2 < Observe that f feRa = if z < o 2 . If 

z > a,2, then the area of the set {/ G R3 : / 2 < ° 2 } is the difference of the area of two triangles (see 
Fig. 3(d)). Direct computations yield 



feR 3 





ai (z — a 2 )(2+a2 — 2) 
k 2(a 2 -l)(z-l) 2 



if z < a%, 
if z > a 2 . 



> 



Let / G Hb. Then 1{< < z} = 1 ^ < ^ /1 > (<ta _ 1 a ) \,_ l) • We have = ^ (a2 _ 1 ° ) 1 (z _ 1) 

If z > g a ^ 1 , then again we compute the difference of the area of two triangles (see 



o-i z < 



02— 1 ' 
Fig. [3(7)} and infer 







ifz< 



a 2 -l ' 

ai(a 2 -l) if «. >, Q2 
2 f3 It Z > — [ , 



where 33 = 1 



(a2 _l)2 (z _ 1) 2- 

By symmetry, the integrals for / G i?4 and / G Rq are obtained by replacing ai with 1 — ai in the 



formulas for R 3 and R 5 (see Fig. 3(c) and 3(g)). Thus, 



f&Re 





(1 — ai)(j8 — a2)(z+Q2— 2) 



2( Q2 -1)( Z -1)2 



if z < a 2 , 

if Z > 02, 





(l-oi)(a 2 -l) 



if Z < 

if z > 



a 2 -l ' 

a 2 
a 2 — 1 ' 



In order to compute ff eR for i = 1,2 we distinguish the cases a 2 < 2 and a 2 > 2. First, let us assume 



that a 2 < 2. Computations in a similar manner as above (area of the shaded quadrilateral in Fig. 3(a) ) 
lead to 



feRi 




if z < a 2 , 

a2) - if Z > 02, 



if z < a,2, 

(1 ' aiK r a2) if^>a 2 - 



<s 



Now let us assume that a 2 > 2. Computing the corresponding area (see Fig. [3(b)] and [3(c)]) yields 



feRi 



feR 2 







if z< 



o a -i ' 
aig 4 if 2 > -^tt, 

if 2 < 




(1 



0,2-1 ' 



where 94 = 1 



(a 2 -l)(z-l)- 



(a) Ri and i?2 
for 02 < 2 



(b) i?x for 
a 2 > 2 



(c) R 2 for 
a 2 > 2 



(d) fl 3 



(e) R 4 



i 



(f) ^5 



(g) Re 



Figure 3: The shaded regions satisfy l{t < z} = 1. 

We can aggregate the formulas for regions R\ and i? 2 , R3 and R4, and R$ and i?6 m order to eliminate 
the parameter 04. It follows 

if z < a 2 and a 2 < 2, 
if z > a 2 and a 2 < 2, 



2 — a2 



/efli /6-R2 



feRi feR2 



if z < -^r and a 2 > 2, 

(a 2 -i)(,-i) if2>^ T anda 2 >2, 



feR 3 feRi 



if z < 0,2, 

(z — a2)(z+a,2—2 
, 2(a 2 -l)(*-l)* 



(z - a2 ) (z+Q2 _ 2 ) j fz>fl2j 



ifz<^-r, 

— a 2 — 1' 



0,2 — 1 



/efl 5 /e-R 6 



(* ~ (a 2 -l)^( Z -l)^) ifz> ^FT 



We are now ready to state our results in terms of the lattice width. Reinterpreting the formulas above 
we obtain for a 2 < 2 (i.e. w = a 2 ) the integrals 

(0 i£z<w, 
IfEf Xz>w, 



feRi feR 2 



feR 3 feRi 



feR 5 feRe 
For a 2 > 2 (i.e. w = -^r) we have 



feRi feR 2 



if z < w, 

(z-w)(z+w-2) jf 



if2!< « . 

It? — 1 ' 

(W-I) 2 (z~lf~l -r W 

. 2(w-l)(z-iy 2 11 z w-i- 



if z < w, 
if z > w, 
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feR 3 feRi 



ilz<-^r. 

— W—X ' 

(w-l) 2 (z-l) 2 -l w 





(z-m)(z+tu-2) 
, 2(u>-1)(>-1) 2 



if Z < 

if z > 



Note that w < — ^ as it; < 2 (see [4] for a proof that any triangle of type 2 satisfies 1 < w < 2). Hence, 
together with ([3]) we arrive at the following formula. 



P T2 (z) > 







Si 



if 1 < 2 < w, 
ii w < z < 



w-l ' 

.91 + .92 if < 2 < 



(7) 



with <?i 



(2— u;)(2'Ujs — lu- 



and ^2 



™ 2 ( 2 -i) a 



Let z > 1 be given. We want to show that the 



probability to improve upon the split closure by a value of more than z when adding a type 2 triangle 
becomes the smaller the closer the type 2 triangle is to a split, i.e. the closer w is to 1. For w being 
sufficiently close to 1 we have w < z < Just substitute 1 for w in ([7]) to infer that P T2 (z) converges 

to 1 for any z > 1. Therefore, 1 — P T2 (z) tends to 0. In other words, the chance that adding a fiat type 2 
triangle improves the split closure by a value of more than z with respect to our strength measure tends 
to 0. In terms of the vertex (01,02) this happens if ai converges cither to 1 (i.e. T2 converges to Si) or 
infinity (i.e. T 2 converges to S%). 

Two special cases are of interest: z — 2 and z = |. Plugging in z = 2 in Q yields 1— P T ' 2 {2) < 4< ^ w ~ i 1 '> 



for any 1 < w < 2 (sec Fig. 4(a) I and plugging in z — | in (JT]) leads to 



5' 





(3-2tu)(4w-3) 



if § < w < 2, 
if 1 < w < I 



(see Fig. |4(b)| . 1 — P T ' 2 (2) is interpretable as the probability that the inequality corresponding to T2 
closes more gap to i?j than the inequality corresponding to any type 1 triangle. On the other hand 
P T2 (|) can be seen as the probability that adding T2 to the split closure is inferior to adding any type 1 
triangle T\ with rays going through the corners of T\ . 



(a) An upper bound for 
1 - P T 2 (2) in terms of w. 




(b) A lower bound for 
P T2 (^) in terms of w. 



Figure 4: Bounds for 1 - P Ta (2) and P T2 (|). 



The interpretation is based on the fact that any type 1 triangle (assuming corner rays) has a strength 
between | and 2. On the one hand type 1 triangles are inferior to flat type 2 triangles by comparing 
the worst case strength (see [?])■ On the other hand a type 1 triangle guarantees a strength of | and 
is therefore superior to a flat type 2 triangle on average. For instance, let T2 be a type 2 triangle with 
w{T 2 ) = 1.1. Then 1 - P T2 (2) < 3.4% and P T2 (|) > 92.5%. Thus, with a probability of less than 3.4%, 
T2 is better than any type 1 triangle T±, and worse with a probability of more than 92.5% if corner rays 
are assumed for T\. We point out that these probabilities are quite close to and 1, respectively, even 
though we used just one split instead of the entire split closure. 
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6 Quadrilaterals 



The vertices of the quadrilateral Q are denoted a = (01,02), 6 = (61,62), c = (01,02), and d = (dx,d 2 ). 
By an affine unimodular transformation, we assume that the point (0,0) (resp. (1,0), (0,1), (1,1)) is in 
the relative interior of the facet with vertices 6 and c (resp. 6 and d, a and c, a and d). We further assume 
< ay < 61 < 1, 1 < 02, 62 < 0, and —62 < 02 — 1. The parameters ax, 0,2, 61, and 62 are assumed 
to be arbitrary but fixed. This implies Cl = - {aa _$_ aiba , c 2 = ^, d l = ^^I^E , 
and d 2 = ^57-ri • One easily verifies c\ < 0, < c% < 1, 1 < d\, < d% < 1, C2 < cfe, and 
Area(Q) = ^(02— 6 2 +di — ci). Under these assumptions it is known that w :— w(Q) = min{a 2 — 6 2 , d% — Cx} 
(see [3]). Thus, without loss of generality we assume w = 02 — 62. We decompose Q into four regions: 
iZi := int(Q n {(si.&a) G M 2 : < x- 2 < ^}), i? 2 := int(Q n {(xi.sa) G K 2 : ^ < a: 2 < 1}), 
i? 3 := int((Q\{ J R 1 U J R 2 })n{(xi,x 2 ) el 2 :0<n< 0}), i? 4 := int((Q \ {R x U i? 2 }) D {(x u x 2 ) G K 2 : 
< x x < 1}), where 6 = ^^^^F^^^ 1 (see Fig. |5J. It is tedious but easy to verify that 

C2 < -^0x ^ d-2 and ai < < &i. As in the case of type 2 triangles, instead of taking the split closure, we 
use a single well-chosen split inequality. 




Figure 5: Decomposition of a quadrilateral. 

In regions Rx and R 2 we use the split Si := {(xx,X2) £ ffi 2 : < 12 < 1} and in regions R3 and R4 
the split S2 := {(£1,2:2) G IR 2 : < xi < 1}. Since the computations of i(Q,Sx) for / G Rx U R% and 
i(Q, S2) for / G i?3 U i?4 are straightforward (see Sections 15.11 - [5751 for an illustration on how to do that 
for type 2 triangles) we only state the results. 



i(Q,Si) 


Location of / 


f2-b 2 

02-/2 
1-/2 

/l-Cl 

1-/1 


/ei?i 
feR 2 
feR 3 
feR 4 



As in Section |5.4[ we now compute the integrals in order to obtain a lower bound for pQ(z). For 
simplicity, let := J* /giJ . l{i < z}df for j = 1, . . . , 4. 

6.1 Regions i?i and i? 2 

Let / G i?i. Then l{f < z} = 1 ^ <z^f 2 > fzj. Observe that J feRi = if ^ < -^f <S> z < 

t(7. If z > w, then we distinguish into > C2 and < C2. If > c 2 , then the area to compute is 
a trapezoid with area Ax ; if -^j < C2 the area to compute is composed of two trapezoids with aggregate 
area A2 , where 

-62 \ / w-(bx- oi) z - 61 ai(z - 1 + 6 2 ) \ 

z-iyv «>-i z ~ l ( a 2 - i)(z - i)y ' 

c W W - (61 - Qi) 01(62 - 1) - («2 - l)6l \ 

/ V w-1 axb 2 - (a 2 - l)6i / 

-62 \ / z fli(6 2 -l)-(a2-l)6i \ 
z — lj\z — 1 ai&2 — (02 — l)6l /' 



Ai = 
A 2 = 



\ 




62 






- 1 


\ 




62 










1 


(1 




f 2 
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We obtain 

if 1 < z < w, 



A 1 iiw<z<z^, 
[A 2 if^-<z<+oo. 



Let / G R». Then l{t < z} = 1 & < z ^ f 2 < S=%. Thus, J /eRa = if < jfe ^ z < w . 

Otherwise, we distinguish into z "^_° 2 < d 2 and z ~_f£ > d 2 . In the first case the area to compute is a 
trapezoid with area A3; in the second case the area is composed of two trapezoids with aggregate area 
At, where 



.4 



3 



1 / ' z — a 2 —b 2 \ / w — (bi — a±) z — 1 + a\ (z — o 2 )(&x — 1) 



2\z-l w-lj\ w-1 z-1 b 2 (z-l) 



A - 1 ( z ~ a ^ d \ ( z , a 2 (l -61) — (1 -01)62 



2\z-l J\z-1 (a 2 -l)(l-6x)-(l-ox)6 2 

1 ( _ -6 2 "\ (vi - (61 - ffli) , 02(1 - bx) - (1 - ax)&2 



We infer 



2 V " w-lj\ w-1 (o 2 -l)(l-6i)-(l-ai)6 2 



if 1 < 2 < K7, 
if W < Z < 




6.2 Regions f? 3 and f? 4 

Applying the same procedure as in the previous section leads to 

If —c\ -ci \ / (a 2 - l)(di - 1) (02 - l)(z — 1 + ci) c 2 c 2 



-4 



5 



2V^x-Ci-l z - 1/ V(l - ai)(di - a - 1) (l-ai)O-l) dx - ci - 1 



1 / -ci \ / (a 2 - 1)(2 - a x ) a x b 2 c x {a 2 - 1) + c 2 (l - a x ) 



2 \di — c\ — 1 / \ 1 — ai &i (1 — ai)(di — ci — 1) 

1 / -ci \ / ai6 2 , O1C2 — cx(aa — 1) 

1 «i 7 02-1 



A 
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2 V z - iy \ °i - 1) 

1 (dx-l)(z-dx + cx) /p2(l-oi) + (a 2 -l)(di-l) , (a 3 -!)(<*! -1) b 2 (z - d x ) 



2 (z-l)(di -ex -1) V (1 - oi)(di - ci - 1) (l-oi)(«-l) bx(z-l) 



1/ -ex ^ ^c 2 (l-ax) + cx(a 2 -l) , (02 -1)(2- 61) -63(1-01) 

A 8 = 7: °1 



2 \ d\ — cx — 1 / \ (1 — ax)(dx — cx — 1) 1 — ax 

1^-dx \ f{a 2 - l)(z-di) 6 2 (dx-l) (a 2 -l)(2-6i)-6 2 (l-ai) 
— ox 



2 V ^ — 1 V V (01 - l)(z - 1) (l-6x)(z-l) 1-ax 
and finally 

if 1 < z < di — Cx , 
5 if dx-cx <z< si^a, 
/ G fl 3 LA 6 if £J^£L < z < +00, 



if 1 < z < di - cx 
A 7 if dx - cx < z < - 
/en, (A 8 if^<z<+c 



6.3 Approximation for P®(z) 

Adding the integrals and dividing by the area of the quadrilateral gives the lower bound for pQ{z) which 
we wanted. Thus, 

p Q ^j^mt j i{*^>d/. ( 8 ) 
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We did not succeed in showing algebraically that this lower bound tends to 1 when w converges to 1. 
However, we performed simulations supporting our conjecture. In these simulations we let w converge 
to 1 for several values for ax, a,2, bx, and 62. Concretely, we discrctized an, a 2 , bx, and 62 within their 
ranges: ax G {0.001, 0.002, . . . , 0.999}, bx G {ax, a 1 + 0.001, . . . , 0.999}, o 2 G {1.999, 1.998, . . . , 1.001}, and 
&2 G { — (ct2 — 1), — (d2 — 1) + 0.001, . . . , —0.001}. In all cases we observed that the lower bound became 
close to one. This is not a proof, but gives strong indication that this should hold in general. 

We point out that under the additional assumption ax = bx (implying a\ = 6 = b\ and C2 = Ej~[ = d 2 ) 
the lower bound ([8|) coincides with the lower bound of a type 2 triangle which is stated in (0 . This can 
be seen by simply substituting ax for 61 in the formulas for Ax, ■ ■ . , Ag. 



COROLLARY 6.1 Let Q be a quadrilateral meeting the assumptions stated at the beginning of this section 
which, in addition, satisfies a\ — bx- Let w := w(Q = ). Then 

{0 if 1< z < w, 

(z — w)(2wZ — W — z) -j? , ^ w /n\ 

Jkz-ir iiw<z<— , (9) 

{z - w)(2wz - w z) + (w -lf( Z -Xf-X i{ ^ <z<+OQ ^ 
w 2 {z—l) 2 w—1 

Moreover, for any z > 1, (z) tends to 1 if w converges to 1. 



7 Type 3 triangles 

By an affine unimodular transformation, we assume that the type 3 triangle T 3 satisfies T 3 n Z 2 = 
{(0,0), (1,0), (0, 1)} and that each facet of T 3 contains one of these points in its relative interior. The 
three vertices are denoted by a = (01,02), b = (61,62), and c = (01,02). We further assume 1 < ax, 
< 02 < 1, < 61 < 1, and 61+62 < 0. Let ax, 02, and 61 be arbitrary but fixed. Thus, the 
other parameters are 6 2 = - "'^Z^ , cj = (ni _ 1)(1 l 1 i°) t> 7i ) a 6 1 1 na(1 _ 6l) , and c 2 = - (ni _ 1)(1 ° 1 °' ) K 1 ) a3(1 _ fcl) ■ 
One easily verifies 6 2 < 0, ci < 0, 1 < c 2 , < ci + c 2 < 1, and Area(T 3 ) = |(ai + 02 — 6 2 — Ci). 
Under these assumptions we have w := w(Ts) = min{c2 — 62,01 — Ci,ai + 02 — (61 + 62)} (see [1]). 
Without loss of generality we assume w — c 2 — 62 < ai — c\ < ai + 02 — (61 + 62). During this 
section we consider the three splits Sx '■= {(^1,^2) G M 2 : < x 2 < 1}, ^2 := {(xi,x 2 ) G K 2 : 
< xx < 1}, and S3 := {(xx,x 2 ) G R 2 : < x\ + x 2 < 1}. We decompose T 3 into six regions: 
Rx := mt(T 3 n {(x u x 2 ) el 2 :0<i 2 < R2 ■= mt(T 3 n {(xx,x 2 ) G R 2 : ^Ej^ < x 2 < 1}), 



R 3 

TO 2 



= int((T 3 \{i?iUi?2})n{(a;i,a ; 2) GR 2 :0<xi < ai I?_ Cl }), Rx := int((T 3 \ {Rx U R 2 }) n {(an , x 2 ) G 



^Eft^ < *x < 1}), ^5 := int((T 3 \ U| =1 ^.) n {(x u x 2 ) G R 2 : < Xl + x 2 < ai+ ~^I&W )' 
and i? 6 := int((T 3 \ U 4 J=l R 3 ) fl {(xi,s a ) G R 2 : ai+ ~^^^ 1 ) +t)3) < »i + ^2 < 1}) (see Fig. We point 
out that i?5 could be empty. It is tedious but easy to verify that 02 < C2 Zx-b 2 < 1» 61 < ai Zii Cl < 1) 
and < -(fri+M < c + c 

ai+a 2 -l-{bi+b 2 ) 1 z 




Figure 6: Decomposition 

For each region, we use a single split inequality to 
i?2 we use the split Sx, in regions i? 3 and Rx the split 



of a type 3 triangle. 

approximate the split closure. In regions Rx and 
5*2, and in regions R$ and Rq the split S3. Thus, 
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in each region Rj we choose a split which covers Rj and the convex hull of (0, 0), (1, 0), and (0, 1). The 
following table states the values for i(T 3 , Si) for the regions R\ to R 6 . 



tiTr, Si) 


Location of / 


C2-/2 
1-/2 

Jl-ci 

01-/1 

/l+A- (61+62) 

/1+/2 

ai+o 2 -(/i+/2) 


feRi 
feR 2 
feR 3 

feR 4 
feR 5 
feR 6 


l-(/l+/ 2 ) 



The computation of the integrals Jf €R . l{i < z}df for j = 1, . . . , 6 gives no new insights as well. We 
only state the results. For simplicity, let ff eR '■= ff<=R 1{* — z }df f° r J : ™ 1> • ■ • 1 6. Let 

1 / — 6 2 —62 \ ( b\ b\ a\ ( c 2 — 1 z — I + bo 

A\ 



2 \W — 1 Z— 1 2—1 1 — (22 \ W — 1 Z — 1 

1 / -6 2 \ /(l -a 2 )6i +ai(c 2 - 1) a 2 6i-ai& 2 ; 
^2 = 77 - a 2 



2\w-\ J \ (1 - a 2 ){w - 1) 

1 / -6 2 \ / 02^1 - (ai - 1)&2 02^1 - (ai + 
77 a 2 



A3 

A, 
A 6 



2 V 2- 1/ V a 2 (z-l) b 2 
1 ( Z — C2 -b 2 \ ( b\ b\{z — c 2 ) ai /cj-l cj-l 



2 \ z — 1 w — 1/ \u> — 1 62(2 — 1) 1 — a 2 \ io — 1 z — 1 
1 / — ci — ci \ / a 2 a 2 (z — 1 + ci) 



2 \ai — c\ — \ z — 1 J \ai —c± — l {ax — l)(z — 1) 
1 a 2 /&i(z-l) + ci 



2 6x(oi-l) V 
1 ( z — a\ — ci \ f a 2 a 2 



2 \ z — 1 ai — ci — 1 / \ ai — ci — 1 z — 1 J ' 
1 / -62 



^4 _ 2 . if ^ 

2 V a i + a 2 - (61 + o 2 ) - 1 

8 2 Vai + a 2 - (61 +62) - 1 62/ V fl i +«2 - (61 +o 2 ) - 1 
. I/61 + 62 61 + 6 2 \ / -6 2 

* = 2 r - * - ~ l^r- 1 

, _1 f 6 2 (z-(ai+a 2 )) / 2 

Al ° - 2 ■ iFWl 

. _ 1 ( ~b\{z - (ai + a 2 )) 6i\ / 6 2 (z - (ai + a 2 )) 
11 ~ 2^ (6i + 6 2 )(z-l) + &JV (&i + &a)(*-l) 

■ _ 1 ai + a 2 - 1 \ / -62 

12 ~2^ 6 2 01 +a 2 - (61 + 6 2 ) - iy Ui +«2 - (61 +6 2 ) - 1 

Al 3 = \-{c2~ I)' , 

A M = |-(£-ci) (C2-1), 

. _ 1 / ai + a 2 - 1 61 \ / -62 1 

15 _ 2 V ai+a 2 -(6i+6 2 )-l 6 2 J \ai + a 2 - (61 + 6 2 ) - 1 

1 ( , v ai + a 2 — 1 \ / 2 — a 2 

Ai 6 = - l-( Cl + C2 )-^_ c 2 -_ f 
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A 1 ~ "2 f-i I I \ a l+ a 2-l 

An = 1 - (ci + c 2 ) 



We obtain 



/ 







if 1 < z < 



w. 



feRi feR 2 



f£R 3 feRi 

'0 




= < 



feR 5 



A 7 - A 8 
A 7 - A 8 



= {Ai+A z \iw<z<2^^, 
A 2 + A 3 if 22-^2 < z <+oo, 



if 1 < z < ai — ci, 
ifai-C! <z<^, 
A 6 if ^f 1 < z < +oo, 

if 1 < z < oi + a 2 - (6i + b 2 ) 
or ai + a2 — 6i — 1 > 0, 
A, if ai + a 2 - (&i + b 2 ) < z < 1 - b 2 
and ai + a 2 — b± — 1 < 0, 

if 1 — b 2 < z < +oo 
and ai + a 2 — foi — 1 < 0, 





^4io 
Am 



An - A12 
An 



Ai 3 -A14 
A13 - A14 



Ai 5 
Ai6 



Ale 
Ai7 



Ai, 



if 1 < z < ai + a 2 — (bi + b 2 ), 
i iai +a 2 -(b 1 + b 2 )<z<^±^±fl 

and ai + a 2 — b\ — 1 < 0, 
i i a 1 +a 2 -(b 1 + b 2 )<z< a -l±^^fl 

and ai + a 2 — b\ — 1 > 0, 

if ai lT ( -! 1+ ^ 2) < z < +oo 

1— (C1+02J — 

and a\ + a 2 — b\ — 1 < 0, 

if ai t a r ( f i+ , C2) < z < +oo 

l-(ci+c 2 ) — 



and ai + a 2 — bi — 1 > 0. 
This leads to the lower bound for P Ta (z) which we wanted, namely 



P Ts (z) > 



1 



6 

Area(T 3 ) £ . 



l{t<2K. 



As for quadrilaterals, we did not succeed in showing algebraically the convergence of this lower bound 
to 1 when w converges to 1. However, simulations with discretized parameters a\ £ {3, 4, . . . , 1 000 000}, 
a 2 E {0.001, 0.002,..., 0.999}, and h G {0.001, 0.002, 0.999} (where b 1 < ai+ ° 2 2 _ 1 ) such that to 
converges to 1 suggest that this is the case. 
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